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We investigate how the nonhnearity of the Zeeman shift for 
strong magnetic fields affects the dynamics of rf field induced 
evaporative cooling in magnetic traps. We demonstrate for 
the ^^Rb and ^^Na F — 2 trapping states with wave packet 
simulations how the cooling stops when the rf field frequency 
goes below a certain limit (for the ^^Rb F — 2 trapping state 
the problem does not appear). We examine the applicability 
of semiclassical models for the strong field case as an extension 
of our previous work [Phys. Rev. A 58, 3983 (1998)]. Our 



recent ^'^Rb 



results verify many of the aspects observed in 
experiment [Phys. Rev. A 60, R1759 (1999)]. 
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I. INTRODUCTION 

Bose-Einstein condensation of alkali atoms in magnetic 
traps was first observed in 1995 [Q, and since then the de- 
velopment in related research has been been very swift. 
Typically the hyperfine state used in the alkali exper- 
iments is the F = 1 state, although condensation has 
been demonstrated for the ^^Rb F = 2 case as well [^. 
The trapping of atoms is based on moderate, spatially in- 
homogeneous magnetic fields, which create a parabolic, 
spin-state dependent potential for spin-polarised atoms, 
as shown in Fig. |l|(a). For slowly moving atoms the trap- 
ping potential depends on the strength of the magnetic 
field B but not on its direction [^. In practice the field 
is dominated by a constant bias field ^bias, which elimi- 
nates the Majorana spin flips at the center of the trap. 

In evaporative cooling the hottest atoms are removed 
from the trap and the remaining ones thermalise by in- 
elastic collisions. This leads to a decrease in temperature 
of the atoms remaining in the trap [0-p[. Continuous 
evaporative cooling requires adjustable separation into 
cold and hot atoms. This is achieved by inducing spin 
flips with an oscillating (radiofrequency) magnetic field, 
which rotates preferably in the plane perpendicular to the 
bias field [g|0] . In the limit of linear (weak) Zeeman ef- 
fect the rf field couples the adjacent magnetic states Mp 
resonantly at the spatial location determined by the field 
frequency [Fig. [l](a)]. Hot atoms oscillating in the trap 
can reach the resonance point and exit the trap after a 
spin flip to a nontrapping state. Using the rotating wave 
approximation we can eliminate the rf field oscillations, 
and obtain the curve crossing description of resonances 
[Fig. Kb)]. 
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FIG. 1. The magnetic trap potentials for spin-polarised 
^^Rb {F — 2). (a) The spin flips that lead to evaporation are 
achieved by an rf field induced multistate transition at a spe- 
cific spatial distance from the trap center. Here Bq — 0.0001 
T. (b) In the curve crossing description the resonances appear 
as degeneracies. Here Bq — 0.0001 T. (c) For strong fields the 
multistate crossing transforms into a sequence on two-state 
crossings between adjacent Mf states. Here Bq — 0.0020 T. 
The circles mark those crossings where the involved adjacent 
Mf states are also coupled. In (b) and (c) we have z^rf = 0.25 
MHz+z/Q. 



The dynamics of atoms as they move past the reso- 
nance point can be described with a simple semiclassical 
model H], which has been shown to agree very well with 
fully quantum wave packet calculations M. The model, 
however, can be applied only if the resonances between 
adjacent Mp states occur at the exactly same distance 
from the trap center. When the nonlinear terms domi- 
nate the Zeeman shifts, the situation changes, as shown 
in Fig. |^(c). The adjacent resonances become separated 
and one expects to treat the evaporation as a sequence 
of independent Landau-Zener crossings as suggested by 
Desruelle et al. in connection with their recent ^^Rb ex- 
periment ||]. We show that there is an intermediate re- 
gion where off-resonant two-photon transitions from the 
Mp = 2 state to the Mp = state, demonstrated in 
Fig. H, play a relevant role. 

In general there is a competition between the adiabatic 
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FIG. 2. The adiabatic potentials (solid lines) and bare Mf 
states (dotted lines) for ^^Na at Bq = 0.0010 T, with rf field 
coupling Q = (27r)20 kHz. The arrow indicates the semia- 
diabatic process for transfer from the Mf = 2 state to the 
Mf = state. Alternatively one can describe the process as 
an off-resonant two-photon transition. 



following of the eigenstates (solid lines in Fig. 0), which 
leads to evaporation, and nonadiabatic transitions which 
force the atoms to stay in the trapping states. In ^^Na 
the nonadiabatic transitions can lead to highly inelastic 
collisions Q. 

In the experiment by Desruelle et al. it was found that 
for a strong bias field the nonlinear Zeeman shifts remove 
some resonances completely, thus making it impossible to 
make a spin flip to a nontrapping state. Our calculations 
confirm this observation. We also show that although 
evaporation could continue via off-resonant multiphoton 
processes, such a process is not practical. The stopping 
of evaporation at some finite temperature occurs for the 
^'^Rb and ^^Na F = 2 trapping states, but not e.g. for 
the ^^Rb F = 2 trapping state. 

In Sec. U we write down the formalism for the Zee- 
man shifts and show the basic properties of the field- 
dependent trapping potentials. We describe the fully 
quantum wave packet approach and corresponding semi- 
classical theories in Sec. Il|, present and discuss the re- 
sults in Sec. |l^, and summarize our work in Sec. ^ 



II. THE ZEEMAN STRUCTURE 

A. 2^Na and ^^Rb 

The Zeeman shifts can not be derived properly in 
the basis of the hyperfine states (labelled by F and 
Mf) poHl^. We need to consider the atom- field Hamil- 
tonian in the (/, J) basis: 



where / and J are the operators for the nuclear and to- 
tal electronic angular momentum, respectively. The first 
term describes the hyperfine coupling; £^hf = ^^hf = 2^4, 
where £^hf is the hyperfine splitting between the F = 1 
and F = 2 states. Here u^ = 1772 MHz for ^^Na and 
UM = 6835 MHz for ^'^Rb. 

The magnetic field dependence arises from the two 
other terms, with C = qjiibB and D = —ajiNB^ where 
the Bohr magneton is iib = eli/2me^ the nuclear mag- 
neton is (In = eh/2mp, and the Lande factor is ^j = 2. 
Here a = 2.218 for ^^Na and a = 2.751 for ^'^Rb. But 
I-^b/ l-^N ^ 1000, and in fact we can omit the third term 
in Eq. (||). 

For 2^a and ^'^Rb we have / = 3/2 and J = 1/2 (lead- 
ing to F = 1 or F = 2 with F = I -\- J). Our state basis 
is formed by the angular momentum states labelled with 
the magnetic quantum number pairs (Mj^Mj). When 
we evaluate the matrix elements of H [using the relation 

/• J = IzJz-\-^{I-\-J--\-I-J-\-)], the states that correspond 
to the same value of Mp = Mj + Mj form subsets of mu- 
tually coupled states. By diagonalising the Hamiltonian 
we obtain its eigenstates. The states which correspond 
to the F = 2 state in the 5 ^ limit (labelled with Mp) 
have the energies Fmf * 

E+2 = ^C, 

E+i = ^\/4A^ + 2AC + C^-A, 

Eo = ^VIa^+cp-a, 



(2) 



E-i 



E. 



-^/aa 



2 - 2AC + C'2-A, 



\C. 



These energies have been normahsed to the energy of the 
F = 2 state for B = 0. In Fig. H(a) and (b) we show the 
Zeeman shifts for all hyperfine ground states of ^^Na and 
^^Rb, but normalised to the ground state energy in the 
absence of hyperfine structure. For small magnetic fields 
(B < 1 T) we get 



^Mf - EmIsMf + (4 ■ 



M]r)e^] 



(3) 



where e = iibB/(2Em). In terms of F and Mp the 
linear Zeeman shift is Em^ = QfI^bBMf = E^sMp as 
the hyperfine Lande factor is gp = 1/2. 

The necessary condition for evaporation is that the rf 
field induces a resonance between the states Mp = 2 and 
Mp = 1. The location of this resonance defines the divi- 
sion between the hot and cold atoms. By decreasing the 
rf field frequency i^rf we both move the resonance point 
closer to the trap center as well as allow more atoms to 
escape the trap. For small B fields all adjacent states are 
resonant at the same location for any u^^. But in case of 
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FIG. 3. The Zeeman shifts for the ground state hyperfine 
states of (a) ^^Na, (b) ^'^Rb and (c) ^^Rb. Note that the situ- 
ations considered in this paper take place in a region located 
very close to B = in the scale of these pictures. 



strong magnetic fields, typically larger than about 0.0002 
T, due to the nonlinear Zeeman shifts the resonances 
separate. Furthermore, the other resonances than the 
Mp = 2 — Mp = 1 one in fact move towards the trap 
center faster, and reach it while the Mp = 2 — Mp — 1 
resonance still corresponds to some finite temperature. 
When Vy-i is lowered further, the other resonances begin 
to disappear. At strong B fields the Mp = state is 
also a trapping state, as shown in Fig. ^, so for effective 
evaporation one really needs to reach the Mp = — 1 state. 
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FIG. 4. The Mf = state becomes a trapping state very 
quickly as B increases. Here we show as an example the situa- 
tion for ^^Na when Bo = 0.0010 T. The oscillation frequency 
for this state is naturally much smaller than for the other 
trapping states. 

At the critical frequency i^cr the crossing between the 
states Mf = — 1 and Mf = disappears. Alternatively, 
for a fixed frequency u^f we have a critical value Bcr for 
the B field; the resonances disappear when B ^ B^r (for 



practical reasons we have chosen to modify B rather than 
i^rf in our wave packet studies). In Fig. ||(a) we show 
the potential configuration when u^^ is slightly below i^cr- 
Since i^cr corresponds to the state separation at the center 
of the trap, it is independent of the trap parameters such 
as the trap frequency. 

For a specific trap configuration i^cr can be converted 
into a minimum kinetic energy required for reaching the 
resonance between the states Mf — 2 and Mf = 1. In 
Fig. |5|(b) we show this minimum kinetic energy in units of 
temperature as a function of magnetic field strength for 
^^Na and ^^Rb, and for the trap configuration used both 
in our simulations and in the experiment by Desruelle et 
al. i. 
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FIG. 5. (a) For B = 0.00155 T and Vri =0.25 MHz+z/q 
the resonance between the Mf — and Mf = — 1 states for 
^^Na can not be achieved (5cr = 0.00152 T). The possibil- 
ity for a tunnelling-like transfer (which could also be called 
off-resonant process via power broadening) exists, though, (b) 
The lowest energy (in temperature units) for which evapora- 
tion is allowed as a function of trap center field Bq for ^^Rb 
and ^^Na with the trap parameters used in this paper. 

In the intermediate region <C 5 ^ ^cr, where the 



necessary crossings exist but are separated, the processes 
take place via two possible routes. We can have off- 
resonant multiphoton processes, that e.g. lead to adia- 
batic transfer from the Mf = 2 state to the Mjr = 
state. This example is demonstrated in Fig. || where we 
show also the eigenstates of the system, i.e., the field- 
dressed potentials. When the relevant resonances are 
well separated, the evaporation takes place via a com- 
plicated sequence of crossings, as indicated in Fig. |l|(c). 
This will be demonstrated with wave packet simulations 
in Sec. Fvl. 



B. s^Rb 

For the isotope ^^Rb we have 7 = 5/2 and J = 1/2, so 
the ground state hyperfine states are F = 2 and F = 3, 
as shown in Fig. ||(c). Now the F = 2 trapping state is 
the lower hyperfine ground state. Thus the behavior of 
the Mp states is different from the ^^Rb and ^^Na case. 
The B field dependence of the states related to F = 2 is 
now 

Q /I 1 . 

E+2 = — - -VdA^nXc + C^, 

Q /I 1 



Eo 



2 2 
3A 1 



V9A^+CP, 



(4) 



'^A 1 / 

E-i = ^- -V9A^-2AC + C^, 



Q /I 1 , 

E-2 = ^- -^9A^-^AC-rC\ 

where now Fhf = 3 A. For ^^Rb we have i^hf = 3036 
MHz. Here the trapping states are now F = 2, M = — 2 
and F = 2,M = -1. If we now define e = (2/3)£ = 
/iB^/(3Fhf) we get approximatively 



Fmp ^ -EmIsMf + (9 - Mi^)e'] 



(5) 



As qf = —1/3, this agrees with the linear expression 
Emf = QfIJ^bBMf = -EmsMf. 

The change of order in the Mf state energy ladder 
means that with increasing B field one never loses the 
crossing points between the adjacent states. In other 
words, if we use an rf field that can couple the states 
Mf = — 2 and Mf = — 1 resonantly at some location xc, 
then we always couple the rest of the states resonantly 
as well at distances larger than xc- In Fig. g we see how 
this leads to a sequence of crossings that allows hot atoms 
to leave the trap without the need for sloshing. One 
must, however, take into account that the kinetic energy 
required to leave the trap is now set by the difference 
between the energy of the Mf = — 2 state at the center 
of the trap, and the energy of the M = (or Mf = — 1) 
state at the point where the states Mf = and Mf = — 1 
are in resonance. In other words, atoms need a kinetic 
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FIG. 6. The effect of nonlinear Zeeman shifts on the evap- 
oration from the ^^Rb F — 2 state. The circles indicate for 
which crossing points the couplings are nonzero. The num- 
bers indicate the order in which an atom moving out of the 
trap traverses the crossings. The kinetic energy required to 
leave the trap is now determined by the difference between 
the trap center and the second crossing. Here Bq — 0.0050 T 
and z/rf =0.25 MHz+z/q, with z/q = (F_2 - E-i)/h. 



energy equal or larger to the energy difference between 
the trap center and the second crossing in Fig. g|. 

In this paper we limit our discussion on the F = 2 case 
only, but it is obvious that for the ^^Rb F = 3 trapping 
states we face the same problem as in the F = 2 case for 
^^Rb and ^^Na. In general for the alkali atoms we can 
expect that the problem will arise whenever we use the 
upper ground state hyperfine state as the trapping state 
at strong B fields. 



C. Trap configuration 

For simplicity we have assumed in our studies the same 
spatially inhomogeneous magnetic field as in the experi- 
ment by Desruelle et al. |9|, except that we have added 
a spatially homogeneous compensation field. This allows 
us to change the general field magnitude (depends on the 
bias field) while keeping the trap shape almost unchanged 
(depends also on the bias field). Thus we set [|| 



B = Bo 



B 



/2 



2B^u 






-^^) + 5^V, (6) 



where B' = 9 T/m, B" / Bub.s = 10^ m-^ and the trap 
center field is defined as Bq = ^bias — ^comp- The actual 
trap is cigar-shaped, which is a typical feature in many 
experiments. We have selected the x direction as the 
basis for our wave packet studies. We set ^bias = 0.0150 
T and use ^comp as a parameter to change Bq. Using 



C = qjjibB with Eqs. (||) and (^) we get the spatiahy 
dependent trapping potentials. 



III. QUANTUM AND SEMICLASSICAL MODELS 

A. Wave packet simulations 

For our wave packet studies we fix the rf field frequency 
to the value v^i = vq +0.25 MHz, where vq = [E-^2{x = 
0) — E-^i{x = 0)]/h. With this setting the atoms need 
typically a kinetic energy about £^kin/^B — 24 /iK in or- 
der to reach the crossing between the states Mf = 2 and 
Mf = 1. With our special definition of u^f the differ- 
ences between ^^Na and ^^Rb appear mainly in the time 
scale of atomic motion (Na atoms are lighter and thus 
move faster), and in scaling of B. For our selected v^f 
we have B^r = 0.00297 T for ^'^Rb and B^r = 0.00152 T 
for ^^Na. We have used the rf field strength Q = (27r)2.0 
kHz), where the rf field induced coupling term is [|U|] 
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in the |F, Mf) basis as indicated. 

The wave packet simulations were performed in the 
same manner as in the previous study [pi. Our initial 
wave packet has a Gaussian shape, with a width of 10 fim. 
For all practical purposes this wave packet is very narrow 
both in position and momentum, and the spreading due 
to its natural dispersion is not an important factor. We 
identify the mean momentum of the wave packet with the 
atomic kinetic energy ^kin, and set £^kin/^B = 30 /iK. 
In the experiment by Desruelle et al. one had typically 
Bq = ^bias = 0.0150 T, which sets the kinetic energy for 
reaching the resonance points (for any practical value of 
z/rf ) too large for realistic numerical simulations. Thus we 
have introduced the compensation field and limit Bq to 
values below 0.0050 T. But the main conclusions from our 
study apply to larger values of Bq and £^kin , and many of 
the results can be scaled to other parameter regions with 
the semiclassical models. 

Another simplification is that we consider only one 
spatial dimension. This is necessary simply because we 
have chosen to work with relatively large energies, such as 
30 /iK. Numerical wave packet calculations at the corre- 
sponding velocities require on the order of 100 000 points 
for both the spatial and temporal dimensions. As dis- 
cussed in Ref. 0, however, this is not a crucial simplifi- 
cation. 

Basically, we solve the five-component Schrodinger 
equation 



ih 






n{x)^{x,t), 



(8) 



The components of the state vector ^(x, t) stand for the 
time dependent probability distributions for each Mf 
state. The off-diagonal part of the Hamiltonian H is 
given Eq. (|^). The diagonal terms are 






MFhUrU 



(9) 



where m is the atomic mass and Umf{x) are the trap 
potentials as in Fig. |l|(a). For states Mf = —2 and 
Mf = — 1 we use absorbing boundaries, and refiecting 
ones for the others. The numerical solution method is the 
split operator method, with the kinetic term evaluated by 
the Crank-Nicholson approach [[ISyi^ 



B. Semiclassical models 

For small magnetic fields the rf field induced reso- 
nances between adjacent states occur at the same po- 
sition, X = XC' In this situation the spin-change proba- 
bility for atoms which traverse the resonance is given by 
the multistate extension H,^ of the two-state Landau- 
Zener model [f6|. We have earlier shown that for the 
evaporation in^^Na F = 2 state at £^kin/^B = 5 /iK and 
small B this model predicts the wave packet results very 
well 0. 

The solution for the multistate problem can be ex- 
pressed with the solutions to the two-state Landau-Zener 
(LZ) model, so we shall begin by discussing the two- 
state case first. Let us consider two potentials, Ui and 
/72, which intersect at x = xc and are coupled by V. 
For strong B, when the crossings are well separated in 
our alkali F = 2 system, V is equal to hft or y^3/2hQ, 
depending which pair of adjacent states is involved [see 
Eq. §]. 

In addition to the coupling F, the relevant factors are 
the speed vc of the wave packet and the slopes of the 
trapping potentials Umf (^) at the crossing. We define 



a = h 



d{U2 - Ui) 



dx 



(10) 



The speed of the wave packet enters the problem as we 
describe the traversing of the crossing with a simple clas- 
sical trajectory, x = vc{t — to) + xq. This allows us 
to enter the purely time-dependent description where 
the population transfer is given by the two-component 
Schrodinger equation 



ih 



*2(t) 



V \( *i(t) 
V avct J I *2(t) 



(11) 



This is the original Landau-Zener theory. In this form it 
is fully quantum and we can obtain an analytic expression 



for state populations Pi and P2 after the crossing. If state 
1 was the initial state, then 



Pi = 1 - exp(-7rA) 
P2 = exp(— ttA) 



A = 



2V^ 
have 



(12) 



Obviously, the Landau- Zener model is only applicable 
when the total energy is higher than the bare-state energy 
at the resonance point. For more details about applying 
LZ theory to wave packet dynamics see Refs. |14,pT 



And now we return to the original multistate problem. 
According to the five-state case of the multistate model 
(see e.g. Ref. Q) the populations Pmf fo^" the untrapped 
states after one traversal of the crossing are 



Pi=4(l-p)p3, 
Po = 6(l-p)V, 

P_2 = (!-#, 



(13) 



where p = exp(— ttA), and A is defined by setting V = 
hQ/2. This assumes that we were intially on state Mp = 
2. We can see that the final population of the initial 
state, P2, is equal to ex.p[—7rhQ'^ / {2avc)] for both the 
two-state and the multistate model if Hamiltonian (0) is 
used. 



IV. RESULTS 

Typical examples of the atomic wave packet evolu- 
tion for the three trapping states are shown in Figs. M 
and g. They demonstrate the sloshing discussed e.g. in 
Refs. [pUmM. The amplitudes of the components de- 
crease as population is partly transferred to another 
state. Similarly new wave packet components can ap- 
pear at crossings. As a wave packet component reaches 
a turning point it sharpens strongly. In Fig. we have 
Po = 0.0018 T, which means that there is no crossing be- 
tween states Mf = and Mj? = — 1. Population transfer 
from the state Mp = 1 to Mf = is weak. The Mp = 
wave packet component has turning points beyond the 
integration space. 

As sloshing continues Stiickelberg oscillations could 
take place as split wave packet components merge again 
at crossings and interfere (for further discussion, see 
Refs. [ p!4|Jl9| ]). However, the wave packet contains sev- 
eral momentum components and thus such oscillations 
are not likely observed, because they are very sensitive 
to phase differences. In our simulations we saw no major 
indication of inteferences. 

In Fig. ^ we track the trap state populations and their 
sum as the wave packet sloshes in the trap and traverses 
several crossings. The magnetic field values are strong 
enough to ensure that the crossings are well separated. 
We can identify when the various crossings take place al- 
though some of them happen simultaneously. The filled 



symbols indicate the corresponding Landau-Zener pre- 
dictions, and we find that the agreement is excellent. 
Some oscillations appear for the ^^Na case [Fig. Q(a)] 
at times between 3.5 ms and 4.5 ms. These may arise 
from Stiickelberg oscillations, but they do not affect the 
final transition probabilities, supporting our assumption 
that in the end such oscillations average out. Note that 
for ^^Na there is no resonance between states Mf = 
and Mf = —1, but for ^^Rb there is and it is seen as 
a stepwise reduction of Ps = P{Mf = 2) + P{Mf = 
1)^P{Mf = 0). 

Near the critical field Per the probability to leave the 
trap via states Mf = —2 and Mf = — 1 varies strongly 
with Pq. When Pq < Per the wave packet meets two 
crossings between the states Mf = — 1 and Mf = as it 
traverses the region around the trap center x = on state 
Mf = 0. At both crossings some population leaks into 
the state M^ = -1, as seen in Fig. |l^ for Pq = 0.0028 T. 
As Pq increases, the two crossing points, on opposite 
sides of X = 0, begin to merge, until they disappear at 
X = when P = Per- Then the transfer between the 
two states becomes off-resonant (or tunnelling), and its 
probability decreases exponentially as a function of some 
ratio of Q and the energy difference between the states 
Mf = — 1 and Mf = at x = 0. This situation cor- 
responds to the parabolic level crossing model |2y]. But 
the main point is that the off-resonant process is unlikely 
to play any major role. 

Finally, in Fig. |ll| we show how the transfer probabil- 
ity between the trap states at the first crossing changes 
as a function of Pq. The multistate process transforms 
smoothly into a two-state process between the states 
Mf = 2 and Mf = 1. The transition zone is rather large, 
though, with Pq ranging from to 0.0010 T. The transfer 
process in this zone is the off-resonant two-photon trans- 
fer demonstrated in Fig. ||. An analogous process can 
occur in atoms interacting with chirped pulses pi]] . 

An interesting point is that the population of the initial 
state is not affected by the fact how the transferred pop- 
ulation is distributed to the other involved states. This 
seems to be typical for the Landau-Zener crossings [ p2[ . 
The solid lines indicate the predictions of the two-state 
model, and the dotted lines the multistate model. They 
change with Pq because the location of the first crossing 
point and thus the wave packet speed vc at this point 
change slightly with Pq. 



V. CONCLUSIONS 

Our results show that in general the semiclassical level 
crossing models offer a clear understanding of the single 
atomic dynamics during the evaporation process. Also, 
we have verified with wave packet calculations that the 
interpretations presented by Desruelle et al. for their 
^^Rb experiment [|[ are correct. The simple picture of 
evaporation at near-zero magnetic fields transforms into a 



complex sequence of two-state crossings at field strengths 
above about 0.0010 T. For all alkali systems where F = 2 
is the upper hyperfine ground state the evaporation will 
stop before condensation as the necessary resonances dis- 
appear too soon as a function of the rf field frequency. 
We have shown that tunnelling does not really play a role 
once the resonances have been lost. Further complica- 
tions arise from the fact that the Mf = state becomes 
a trapping state. 

In experiments, as suggested by Desruelle et al., one 
could avoid the problem by coupling the F = 2, Mp = 2 
trapping state to the F = 1, Mp = 1 nontrapping state, 
or by using several rf fields of different frequencies within 
the F = 2 hyperfine manifold. Although for ^^Rb one 
has observed a long-lasting coexistence of F = 1 and 
F = 2 condensates, theoretical studies [^ predict this 
difficult for ^^Na due to destructive collisions. Thus the 
first approach may apply better for ^^Rb than for ^^Na. 

We have calculated earlier g] that for ^^Na the colli- 
sions between atoms in the Mp = and Mp = 2 states 
are very destructive, with a rate coefficient on the order 
of 10~^^ cm^/s. For practical bias field strengths the Mp 
state is also a trapping state. Thus the efficiency of evap- 
oration is reduced, and the time the atoms spend on the 
Mf = state increase, making it more likely to have a de- 
structive, energy releasing collision. So far condensation 
on the F = 2 state for Na has not been achieved. Even 
in the weak B field case evaporation can produce atoms 
on Mf = state via nonadiabatic transitions. Thus the 
role of inelastic collisions is expected to be enhanced for 
the field strengths considered here. 

Once condensation is reached, however, the nonlinear- 
ity of the Zeeman shifts can be an asset rather than a 
nuisance. For instance, one could create a new type of 
binary condensates by making a selective transfer of part 
of the condensate from the F = 2, Mf = 2 state to the 
F = 2, Mf = 1 state, either by using resonant or chirped 
rf field pulses. Alternatively, two rf pulses of different 
frequencies or perhaps a single chirped pulse might allow 
one to transfer the condenstate from the F = 2, Mf = 2 
state to the F = 2, Mf = state and let it expand nor- 
mally, without the need to switch the magnetic fields off. 
Of course, this would work only when B is so small that 
the trapping nature of the Mf state is not too strong. 
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FIG. 7. The evolution of the wave packet components 
|^Mp(x,t)|2 for ^^Na at Bo = 0.0018 T and Q = (27r)2 kHz. 
(a) Mf = 2, (b) Mf = 1, and (c) Mf = 0. 
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FIG. 8. The evolution of the wave packet components 
|^Mp(x,t)|^ for ^^Rb at Bq = 0.0028 T and Q = (27r)2 kHz. 



(a) Mf = 2, (b) Mf = 1, and (c) Mf = 0. 
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FIG. 9. The time evolution of the trap state populations. 
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FIG. 10. The time evolution of the trapped population for 
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FIG. 11. The population transfer after the first crossing for 
^^Rb as a function of Bq. The solid lines are the semiclassical 
two-state predictions, and the dotted ones are the multistate 
model predictions. The symbols represent the wave packet 
results. 
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